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Abstract 

We derive explicit expressions for the sum rules of the eigenvalues of inhomoge- 
neous strings with arbitrary density and with different boundary conditions. We 
show that the sum rule of order N may be obtained in terms of a diagrammatic 
expansion, with (iV — l)!/2 independent diagrams. These sum rules are used 
to derive upper and lower bounds to the energy of the fundamental mode of 
an inhomogeneous string; we also show that it is possible to improve these ap- 
proximations taking into account the asymptotic behaviour of the spectrum and 
applying the Shanks transformation to the sequence of approximations obtained 
to the different orders. We discuss three applications of these results. 

Keywords: Helmholtz equation; inhomogeneous string; perturbation theory; 
collocation method 



1. Introduction 

In this paper we consider the problem of obtaining exact results for the 
spectral zeta functions of inhomogeneous strings at positive integer values: 

^(^) = El?7 ' ^ = 1,2,... (1) 

n=l " 

where En are the eigenvalues of the string, obeying different boundary conditions 
at the extremities. Remarkably, the exact calculation of Z{s) for integer values 
of s does not require the explicit knowledge of the spectrum, as we have recently 
pointed out in Ref.[H. However, the main focus of Ref.Q was on the possibility 
of obtaining a perturbative expansion for Z(s) for s > d/2 (d is the number 
of dimensions) and use it to evaluate the Casimir energy of the system via 
its analytic continuation. Here we restrict our analysis to the one dimensional 
problem, aiming at obtaining an explicit expression for the sum rule of arbitrary 
integer order n. The extension of these results to higher dimensions is discussed 
in a companion paper Q. 
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As we have pointed out in [l| there is a good number of examples in the 
literature where sum rules have been obtained for different problems H, 3, H , |1, 
H i, S ES El, El, E, [3 [II El, 113, [II : in particular Berry [| and Crandalltial 
have used the sum rules for the modes of certain two dimensional billiards and 
of one dimensional quantum problems respectively to obtain approximations to 
the fundamental mode of these systems. For the case of a quantum bouncer, 
Crandall has obtained the energy of the fundamental mode with a fractional 
error of 10~^, using the sum rule of order 30. 

In this paper we show that it is possible to obtain an explicit expression 
for the sum rule of a string with arbitrary density, and provide a simple set of 
diagrammatic rules which allow one to write down the expression. The actual 
calculation can be conveniently carried out with the help of a computer and 
therefore sum rules of high order can be efficiently calculated. We show that 
very accurate approximations for the energy of the fundamental modes can be 
obtained using a sequence of sum rules. We discuss the application of these 
results in three examples. 

The paper is organized as follows: in Section [5] we derive the general ex- 
pression for the sum rule of order n and state a set of diagrammatic rules for 
Z{n); in Section [3] we discuss the application of the general results of Section [5] 
to three non-trivial problems; finally, in Section 2] we state the conclusions. 



2. Exact sum rules 

In this section we derive exact sum rules for inhomogeneous strings subject 
to different boundary conditions. We consider strings of length a {\x\ < a/2) 
and density T,{x) (S(a;) > 0). 

The Helmholtz equation for the inhomogeneous string is 

(f 

- ^^"(^) = En'S{x)ip„{x) , (2) 

where En and tpn{x) are the eigenvalues and eigenfunctions of the string. At 
the extremities of the string, x = ±a/2, the eigenfunctions tpn{x) fulfill Dirichlet- 
Dirichlet, Neumann-Neumann, Dirichlet-Neumann, Neumann-Dirichlet or periodic- 
periodic boundary conditions. 

Defining (t)n{x) = ^/^{x) ipn{x), as discussed in Ref. eq.([2]) becomes 

\ 1 

-— ^ —==(j)nix) = En(l)„{x) . (3) 



Therefore En and (f>n{x) are the eigenvalues and eigenfunctions of the her- 
mitian operator O = , ^ I — Tr-n- I , : the inverse operator is O^^ = 



2 



Alternatively we may define S,n{x) = (^^53) ~~r' — 4'n{x), which fulfills 
the equation 



In this case En and ^n(^) are respectively the eigenvalues and eigenfunctions 
of the hermitian operator Q = {— J f — j ; its inverse operator 

Observe that 

• the operators O and Q are isospectral and solving any of the eqs.©,® 
or is equivalent to solving the remaining two equations; 

• the eigenvalues of O and Q grow quadratically in n for 

• the spectrum of the inverse operators is bounded both from above and 
from below, < -g- < for n = 1, . . . , 00; 

• the trace of a hermitian operator is invariant under unitary transforma- 
tions and therefore 

n " 

for s = 1,2,...; because of this invariance the trace may be evaluated in 
any basis, in particular in the basis of an homogeneous string, |n): 

Z{s)^Y.^n\[d-'\\n) (6) 



• because of the asymptotic behavior of En for n — > 00, the trace above is 
finite for s = 1, 2, . . .; 

For practical purposes it is convenient to express Eq. ^ directly in terms of 
the Green's functions of the negative ID Laplacian, subject to different bound- 
ary conditions. 

The one dimensional Green's functions have the general form 

G{x,y)^G-{x,y)6{y~x) + G+{x,y)e{x~y) , (7) 

where 

G+{x,y) = G-{-x,-y) , G+{y,x) = G-{x,y) . 

In [Appendix A| we report the explicit expressions for the Green's functions 
for different boundary conditions. 
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Using these Green's functions we are able to write the spectral zeta function 



at s = 1, 2, . . . as: 



/a/2 f-ojl rOj/2 /*a2 

dxi / dx2 ■ ■ ■ dxn-1 / dx. 

-a/2 J-a/2 J -a/2 J -a/2 



3/2 J-a/2 J-a/2 J-a/2 

■ G{XI,X2)G{X2,X3) . . .G{x„-i,Xn)G{Xn,Xi) S (.Xi ) . . . S (x„ ) . (8) 

An equivalent expression for Z{n) is easily obtained using the "x-ordered" 
product of the Green's functions: 



/a/2 /■xi rXn-2 rx„-i 

dXi / dX2 ■ ■■ dXn-1 / dXn 

-a/2 J-a/2 J-a/2 J-a/2 

a(a;i,...,x„) S(xi)...S(a;„) , (9) 



where we have defined the n-point Green's function 

g{xi, . . . ,Xn) = [G{X1,X2)G{X2,X3) . . . G{Xn-l , Xn)G{Xn , Xi)]j, (10) 

and 



[f{xi,X2,...,Xn)]-p= f^^- 
permutations 



PU • • • ) -^Pn 



For example: 

[/(a;i,X2)]p = f{xi,X2) + f{x2,Xi) 
[f{Xl,X2,X3)]p = f{xi,X2,X3)+f{xi,X3,X2) + f{X2,Xi,X3) 
+ fix2,X3,Xi)+f{x3,Xi,X2) + f{x3,X2,Xi) . 

We write explicitly G up to order 5: 
g{xi) = G+{xi,xi) 

G{X1,X2) = 2 [G+{xi,X2)f 
gixi,X2,X3) = 6 G+{xi,X2)G+{xi,X3)G+{x2,X3) 

g{xi,X2,X3,X4) = 8 [G+(xi,a;2)G+(a;i,a;4)G+(a;2,X3)G+(x3,a;4) 

+ G+ {Xi ,X3)G+ (Xi ,X4)G+ {X2 ,X3)G+ {X2 , Xi) 

+ G+ {xi ,X2)G+ {xi , X3)G+ {x2 ,Xi)G+ (xs, xa)] 

g{xi,X2,X3,X4,X5) = 10 [G+(a;i,X4)G+(.Ti,.X5)G+(a:2,a;3)G+(x2,X5)G+(:r3,a;4) 

+ G+(xi,X3)G+(.'i;i, .'i;5)G+(x2,X4)G+(a;2, a;5)G+(a;3, X4) 

+ G+(xi,X2)G+{xi,xr^)G+ {x2 , X4 ) G+ (.X3 , .X5 ) G+ {xs , X4 ) 

+ G+(a;i,a;2)G+(a;i,X4)G+(x'2,x'5)G+(x3,x-5)G+(x3,x4) 

+ G+ {Xi ,X2)G+ (Xi , X5)G+ {X2 ,X3)G+ {X4 ,X5)G+ (xs , X4) 
+ G+(X1,X2)G+(X1,X3)G+(X2,X5)G+(X4,X5)G+(X3,X4) 
+ G+(xi,X4)G+(Xi,X5)G+(X2,X3)G+(X2,X4)G+(X3,X5) 
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Figure 1: Diagrams for Gi^i), 6(2^1 , 2^2) and G{xi, X2, X'^). 

+ G+ (Xi , X3)G+ {Xi , X4)G+ {X2 ,X4)G+ {X2 ,X5)G+ (X3 , 0:5) 

+ G+ {xi , X3)G+ {xi ,X5)G+ {x2 , a;3)G+ (2:2, Xi)G+ [x^ , 2:5) 

+ G+(xi,x^)G+{xi,Xi)G+{x2,x:i)G+{x2,x^)G+{xi,xz) 

+ Gj^{xi,x2)G+{xi,XA)G+{x2,x^)G+{x'i, x^)G+{xa,x^) 

+ G+{xi,x2)G+{xi,X'i)G+(x2,XA)G+{x'i, xz)G+{xi, x^)] 

We can therefore calculate Z{n) with the diagrammatic rules: 

• Draw n points xi, . . . , x„ on a line; 

• Connect each point to any two other points in all possible inequivalent 
ways excluding the disconnected diagrams and the diagrams corresponding 
to a cyclic permutation of the points; 

• Associate a density S(xi) at each point Xi (i = 1, . . . , n); 

• Associate a factor G+{xi^Xj) to each line connecting Xi to Xj {i < j); 

• Multiply the result by a factor 2n, corresponding to the n cyclic permu- 
tations of each inequivalent configuration and to the 2 possible directions 
in which each diagram can be traveled; 

• Integrate the expression obtained from the steps above over the internal 
points: 




It is easy to convince oneself that working to order n there are n!/2n = 
{n — l)!/2 independent diagrams (for n > 2). For example in FiglT]we plot the 
diagrams for Gixi), Gixi,X2) and G{xi,X2,X3); in Figl2]we plot the diagrams 
for G{xi,X2,X3,X4): in this case there are 4!/8 = 3 inequivalent diagrams. 

We can now use any of the Eqs.®, (|S]) or ([5]) to evaluate the sum rules for 
the string, although Eq.® is to be preferred for practical purposes, because of 
its simple diagrammatic representation 0. 



^In Ref. m we have used Eq.(j6l to evaluate explicitly the sum rules of certain strings. 
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Figure 2: Diagrams for Q{xi,X2,X3,X4). 



It is particularly simple to evaluate Z{1), which reads 

Z<™>(1) = (^± + !^'^EMdx (12) 

+a/2 



= / " (13) 

= r'\l-x)n^)dx (14) 



a/2 
+a/2 

-a/2 

z(^^)(i) = r"^y:{x)dx. (15) 

Therefore we obtain the interesting relations: 

ZiDD+NN)^^^ = Z(^^)(l) + Z(^^)(l) = y(E) (16) 

ZiDN+ND)^^^ = Z(^^)(l) + Z(^^)(l)=a2(i]) (17) 

= ^(S), (18) 

where (S) = ^ J^^^j^ T,{x)dx is the average density of the string. Apparently 
these simple relations have not been found previously. 

The diagrammatic rules that we have enunciated allow one to obtain the 
spectral zeta functions of a string at large integer values using a computer; as 
we shall see soon, these sum rules are useful since they allow one to derive 
approximations to the energy of the ground state. 

As early as 1776, Waring in his " Meditationes Analyticae" considered the 
roots a, /3, 7, . . . of an equation (with a > /3 > 7 > . . .) and discussed the 
convergence of (a^" + + 7^" + . . towards the largest root 00. This 



■^This is originally reported by Whittaker and Robinson in Ref. |19|| and by Watson in 
Ref. 1 2011. This result appears at pag. 375 in the edition of 1785, freely available at google 
books I21I. 
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result may be expressed as 0] 



El = lini Z(s)-i/^ . 

S— >-oo 

Little later, in 1781, Euler used the inequalities 

Zis) 



Z{s)-^/' <Ei< 



Z{s + l) 



(19) 



(20) 



to approximate the lowest Bessel zero (see also [HIMIl)- 

Berry [| and Crandall ^ have used these formulas to obtain precise ap- 
proximations to the energy of the ground state of certain systems: in particular, 
Berry in Ref. jsj has introduced a " semiclassical zeta approximation" , which 
amounts to evaluate the energy of the ground state with the formula 



where, following Berry's notation. 



N 



Z{2) 



oo 



/ . p2 
j=N+l 3 



Z^{2) + Z+{2) . 



(21) 



(22) 



Here Z^{2) is estimated using a semiclassical approximation for the higher part 
of spectrum. In this way Berry has obtained the lowest energy of certain two 
dimensional billiards within a 1% error. 

The approach followed by Berry is particularly convenient in our case, both 
because we have at our disposal explicit formulas for the spectral zeta func- 
tions at integer positive values, and because we know the first few terms of the 
asymptotic behavior of the spectrum of the string 

We may generalize Eq.([2T|) to 



El = [Ziq) - Z+{q)] 



-1/9 



(23) 



where q = 1,2, . . notice that since Z{q) > Z^{q), one recovers Eq. (fT9|) in the 
limit q oo, neglecting Zf{q). 

The procedure followed by Berry may also be applied to the calculation of 
the eigenvalues of the excited states: in this case the energy of the -n}^ state is 
obtained as 



En = 



Z{q-n + l)-Z+{q-n + l)-Y, 



1 



E'^ 



-n+l 



-l/(g-ri+l) 



(24) 



•^The asymptotic behavior of the spectrum of an inhomogeneous string with Dirichlet 
boundary conditions is En ~ Axn^ + A2 + A-j/n^ -|- . . .. The coefficients Ai and A2 may 
be obtained using the WKB method; in Ref. |22| we have discussed a method which allows one 
to obtain A3 and higher order coefficients. 
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For example, the second eigenvalue is obtained using Eg.ip^ and it reads: 

E2=\z{q-l)-Z+{q-l)~{Z{q)-Zt{q)t-'^'''y"'^'''^ • (25) 



Notice that as long as the Z^{q~n+\) are evaluated exactly, these equations 
are also exact; for sufficiently large q (with q 3> n), Z^{q — n + 1) may be 
estimated accurately using the asymptotic behavior of the spectrum. 

In |Appendix B| we show that, for n — > oo, for any of the boundary conditions 
considered 



En -> Oien + P 

_ , J-a/2 16S(a;)V2 



c-a/2 4S(x)T." (x)-bT.' (x f , 
J-a/2 16S(a;)V2 



(26) 



where e„ are the eigenvalues of a homogeneous string of length a and a{x) = 
\/S(2/)(iy. Notice that, for the case of Dirichlet-Dirichlet boundary con- 
ditions, this formula provides the coefficients Ai and A2 of ref422]. 

Thus we are interested in calculating the semiclassical approximation to 
Z+{s), which we can Z+{s): 

00 

Z^is) = E ("^J- (27) 

j=n+l 



For instance, for s = 1: 



, , acothf^) T 

^(..H.. _ «tanh(|^) «coth(|^) , 
atanhf^i^') 



a coth 



1 



(32) 

Explicit expressions for s = 2,3, .. . can be also obtained, although we do 
not report them here. 

When we use the formulas above setting /? = we obtain the interesting 
relations: 

3a 3 
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ziDN+ND)+^^^ = '^ = a{a/2f (34) 
a 

(1) - 12a- 12 ^^^^ 

which have the same form of eqs. dTSl) . ((T7)) . apart from (T(a/2)^ ^ a^(E). 

We wiU now discuss a different method to obtain accurate approximations 
to the energy of the fundamental mode of an inhomogeneous string of arbitrary 
density. As we have seen, it is possible to obtain explicit expressions for the sum 
rules Z(ri) corresponding to different boundary conditions: the only limitation 
to the calculation of Z{n) is the factorial growth of the number of terms in Q 
contributing to Z{n) when n is large and the difficulty in performing analytically 
the integrations over the coordinates xi, . . . , a;„. 

If we assume that this complications may be overcome up to some N , then 
one has a sequence of sum rules, Z{1), Z{2), . . . , Z{N), which can be used to 
obtain a sequence of approximations to Ei, as explained before. 

It is easy to convince oneself that the terms in this sequence converge expo- 
nentially to Ei: for example, we notice that for s 3> 1, we have 

^l/s - El f El 



« i?i - (^-^ j +... (36) 

where we have only kept the leading correction. A similar behavior can also be 
inferred for the sequence of ratios Z{s)/Z{s + 1): 

Z{s) f EiY E2-E1 

Notice that the lower bound is more accurate than the upper bound. 

Sequences with transient behavior of this kind can be efficiently extrapolated 
using the Shanks transformation [i^ : in this way, starting with the sequence of 
values of Z(n)~^/" one obtains a new sequence 

Z{n - l)-^Z{n + - Z{n)-i 

Z{n - + Z{n + 1)"^ - 2Z{n)~i 

which converges more rapidly to Ei. The new sequence has N — 2 terms. 

Notice that the Shanks transformation can be applied repeatedly, as long 
as the sequence at one's disposal has at least three terms. Therefore, one can 
in principle obtain a large gain in precision by eliminating several transient be- 
haviors from the original sequence. The advantage of this procedure is that we 
are dealing with exact sum rules and therefore we do not have to worry about 
round-off errors which would necessarily be present in a numerical calculation: 
moreover, the Z(n) are calculated explicitly as functions of the physical param- 
eters in the problem and therefore the Shanks transformation will also provide 
an analytical expression in terms of the physical parameters. 

We will see several applications of this method in the following section. 
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3. Applications 



3.1. Isospectral strings 

Isospectral strings are strings with different densities, but with the same 
spectrum. A well known example was discovered long time ago by Borg (2^ : 
this string has a density 

with a > —1. Borg proved that for a > —1, all the strings have the same 
Dirichlet spectrum of a string of constant density, corresponding to a = 0. 

Using the sum rules obtained before it is easy to see that these strings 
are only isospectral to the uniform string for Dirichlet boundary conditions. 
As a matter of fact already for s = 1, the sum rule for Neumann- Neumann, 
Neumann-Dirichlet, Dirichlet-Neumann and periodic-periodic all depend on a: 

. £±i (40) 

= i(2a + 3) (41) 

Z(^^)(l) = "(" + ^^ + ^ (42) 
^ ^ 36(a + l) ^ ' 

This result can be better understood noticing that the average density of 
the string depends on a: (S) = p^f^- 

Gottlieb has proved in Ref.[23| that, given a string of length a and with 
density S(a;), the strings with density 

S(x)=C'(2^)'S(e(a;)), (43) 

with 

, aa(a + 2x)+Ax , ^ ^. 

fix) — — ^ (44) 

are isospectral to the first string for Dirichlet be. Notice that ^(x) maps the 
interval (—a/2, a/2) onto itself. The case discussed by Borg is a special case of 
ea. (P5|) and corresponds io T.{x) = 1. 

It is now easy to check the isospectrality of the strings with Dirichlet be: it 
is essential to notice that for the transformation of eq. p4)) 

G^^''\x,y) = Gf ^^(C(x),C(2/)) ^i'{x)^'{y) . (45) 
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Using this property we may express the spectral zeta function for the string 
with density ea.(|43p at arbitrary integer values n as 



/a/2 i-xi (■a;„_2 px„-i 

dxi / dx2 . . . / dXn-1 / dXn 

-a/2 J--a/2 J-a/2 J -a/2 

a/2 nil r£.n-2 ri-n-i 

dxi / • • • / rf^j-i / dS,„ 

-a/2 J-a/2 J-a/2 J-a/2 

■ a(^^)(Cl,.-.,en)S(Cl)...S(en)=4''''^H , (46) 

where we have used the notation Z^,^^"* (n) and Z^^'^-' (n) for the spectral zeta 

functions of the strings with Dirichlet be and with density S and S, respectively. 

This results holds for arbitrary integer n and arbitrary real a > —1 and it 
is consistent with the isospectrality of the two strings. 

It is straightforward to see that the Green's functions corresponding to dif- 
ferent boundary conditions do not obey the transformation (|45l) and therefore 
they are not isospectral. 

Notice that starting from order 2 the sum rules for Neumann be contain a 
non polynomial dependence on a. For example ^^^^^(2) is 

^{NN} = J. liQ^s _,_ 12^7 + 93q,6 _,_ i422a^ + 60210^ 

810a^{a + 1)2 L 

+ 12420a^ + 14220a2 + 8640a + 2160] 

2(a + l)(a + 2)(a(a + 2) + 2) 



3a5 



log(a + 1) 



We do not report here higher order sum rules because of their lengthy expres- 
sions. 

In Fig [3] we plot the bounds for the energy of the fundamental mode of the 
Borg string with Neumann boundary conditions as function of the parameter 
a. The shaded area is the allowed region. The bounds are obtained using the 
sum rules of order 3 and 4. 



3.2. An exactly solvable string 
We consider a string of density 



which belongs to a family of inhomogeneous strings first studied by Horgan 
and Chan ,2^]. The frequencies of these strings can be calculated with arbitrary 
precision, since they are solutions to a trascendental equation. 



In Ref. [24I we have calculated the first 10000 Dirichlet eigenvalues of the 



string (|47)). each with a precision of 200 digits. Using these numerical results, 
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Figure 3: Bounds for the energy of the fundamental mode of the Borg string with NN boundary 
conditions as function of a. 



we have extracted the leading asymptotic behavior of the Dirichlet spectrum of 
this string 

_ 2 2 3 165 73179 81997443 

n^TT n + - - 5^27r2n2 ^ 819207r4n4 ~ 146800647r6n6 ^ ' 

The first three coefficients of this expansion were also obtained analytically 
using a WKB-perturbation expansion, developed in Ref. [2§| . 

We can now use these numerical and analytical results to test our sum rules; 
using Eq. ^ with the help of Mathematica we obtain the exact sum rules for 
Dirichlet boundary conditions: 

Z(^^)(2) = -| + llog^(2) 

^(DD)(3) ^ A iog3(2) + l log(2) 

Z(^^)(4) = + 1^ log4(2) - - log2(2) + — log(2) 

^ ' 46080 81 ^ ^ ' 27 ^ ^ ^ 864 ' 

(^^) ^ _ log5(2) + ^ log3(2) 

^ ' 589824 243 ^ ^ 162 ^ ^ ^ 
475 , 2/~s 917 



log^(2) - — — log(2) 



5184 ° ' ' 27648 

„inm, ^ 11466667 64 , 14 , 4, , 95 , , 

Z(DD>(6) = + log^(2 og^(2 + log^(2 

^ ' 2752512000 729 ^ ^ ' 81 ^ ^ 1296 ^ ^ ' 
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1897 , 2/„N 13183 , 



34560 ° ' ' 368640 

^ ' 31708938240 2187 ^ ^ ^ 729 ^ ^ ' 11664 ^ 

6713 , 3,„, 463043 , 2,^^ 728683 



log (2) + nn.o.on (2) " ,.^..annn log(2) 



103680 ° ^ ' ' 9953280 ° ' ' 147456000 
^ _ 448469829001 + ^ _ jgj 

^ ^ 466121392128000 6561 ^ ^ ' 2187 ^ ' 

95 , 3857 , 4,„, 92749 , 

+ log (2) + log (2) log (2 

2187 ^ ^ 58320 ^ ^ ' 1866240 ^ ^ ^ 

8508391 , 2.„^ 8136221 , 

+ log (2) H log(2) 

5225472000 ^ ^ 1486356480 

„tnm, ^ 5652867433 512 , n, , 56 , 7, , 95 , , 

Z(DD>(9) = log^(2) + log^(2 log^ 2 

^ ' 60881161420800 19683 ^ ^ ' 729 ^ ^ ' 2916 ^ ^ ' 

4837 , 13261 , 4,^, 7288553 , 

- log 2 H log^ 2)H log^(2 

77760 ^ ^ ' 276480 ^ ^ ^ 2322432000 ^ ^ ^ 

65171959 , 2.„^ 1880324004961 , 

- log 2)H log 2) . 

5945425920 ^ ^ ' 699182088192000 ^ 

The factor log(2) in these expressions is related to the average density of the 
string, (E) = J^^^^^^I]{x)dx — |log(2). Notice that the largest sum rule obtained 
here corresponds to a diagrammatic expansion involving 20160 inequivalent di- 
agrams. 

In Table [1] we report the error Z''^^^{q) — zium\q), where 

^LmHl) = X! (num)9 ^ X! p(asym)9 ^^^^ 
fc=l ^fc fc=10'i + l 

and are the numerical Dirichlet eigenvalues previously calculated in 



Ref. [22I and i;'''"^") are given by Eq.®. 

In Table [2] we report the estimates of e[^^^ using Egs. pJ)) and p3|. where 
Z^^^^\q) is approximated with Z^^'^^\q). The underlined digits are exact. 

In Table [3] we report the estimates of ijj;'^'^-' using repeated Shanks trans- 
formations of the sequence in the second column of Table [2] corresponding to 
Eq.dlll). 

In Table |4] we report the estimates of E^^^^ using repeated Shanks transfor- 
mations of the sequence in the third column of Table [2] corresponding to Eg . (|23|) . 
The energy of the fundamental mode is obtained with 17 digits of precision. 

3.3. A string with rapidly oscillating density 
We consider a string with density 

E(x) = 2 + sin(^-(- + ^/^)) , (50) 
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Table 1: - z'Sf '(g). 



q 




^num (ti) 


1 


6.8 




2 


1.2 


10-58 


3 


1.5 


10-67 


4 


1.8 


10-76 


5 


2.1 


10-85 


6 


2.3 


10-94 


7 


2.5 • 


10-103 


8 


2.6 • 


10-112 


9 


2.8 • 


10-121 



Table 2: Estimates of e[ ' using Eq.(lT9]l and Ea.|[23t. 



q 




(z(--)(,)-zr 


1 


6.13866459 


10.22002206 


2 


9.80124983 


10.21851148 


3 


10.15503866 


10.21820809 


4 


10.20660399 


10.21813692 


5 


10.21580556 


10.21811931 


6 


10.21762510 


10.21811486 


7 


10.21800650 


10.21811373 


8 


10.21808942 


10.21811344 


9 


10.21810790 


10.21811337 
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Table 3: Estimates of i?J using repeated Shanks transformations of the sequence in the 
second column of Table [2] corresponding to Eq. II19I I . 



Si 


^2 


^3 


^4 


10.19286707426 
10.21540206670 
10.21780418009 
10.21807358764 
10.21810765046 
10.21811245123 
10.21811319374 


10.21809078335 
10.21810761972 
10.21811258058 
10.21811323885 
10.21811332959 


10.21811465291 
10.21811333956 
10.21811334410 


10.21811334408 



Table 4: Estimates of i?J ' using repeated Shanks transformations of the sequence in the 
third column of Table [2] corresponding to Eq. I|23|l . 



^1 


S2 


^3 


^4 


10.2181318565099322 
10.2181151161288641 


10.2181133603626791 
10.2181133457026718 


10.2181133446673210 
10.2181133446659633 


10.2181133446659408 


10.2181135270306264 
10.2181133642743735 
10.2181133468298224 
10.2181133449084580 
10.2181133446933714 


10.2181133447356183 
10.2181133446706434 
10.2181133446662585 


10.2181133446659411 
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with e — )■ 0+ and < 1/2. This example was studied by Castro and Zuazua in 



Ref . [27| , obtaining the asymptotic behavior of the Dirichlet spectrum for e — >■ 0+ 
usin g th e WKB method. More recently, these results have been reproduced in 
Ref. [22|, using an alternative approach, developed by the author: in particular, 
the energy of the fundamental mode, which had been calculated by Castro and 
Zuazua to order e'^ has been obtained in Ref. [22[ to order and reads 



E^^ « ^-^.^ + i.sin^ne^-l^e^ 



64 4 \eJ 1024 

7r(5sin(^) - 116cos(^) +116) , .fin 

In Fig|4]we plot the upper and lower bounds obtained with the sum rules of 
order 4 and 5 for the energy of the fundamental mode of the string ([50]) with 
Dirichlet boundary conditions as function of e. In Figl5]we compare the exact 
asymptotic behavior of e[^^^ for e ^ 1 of Eq. ((5T|) with the approximation 
obtained using the Shanks transformation 

2-(r'r')(3)-i/32-(i5i5)(5)-i/5 _ ^{00)^^-^-1/2 

S — 



Z(-D£')(3)-i/3 + Z(^'D)(5)-i/5 _ 2Z(-D^)(4)-i/4 

In particular for e — 0^ we have 

S « 4.9347 -0.15436^ 

-3 ^ 0.3929- 0.3929 cos ^^^^ - 0.1463e'' 

0.0155 sin - 0.3605 cos + 0.3605 ) (52) 

which should be compared with the exact asymptotic formula 
e[""^ « 4.9348 -0.15426^ 

6^ ( 0.3927 - 0.3927 cos ( — I ) - 0.14466" 



+ ^0.0153 sin - 0.3559 cos + 0.3559^ + . . . (53) 

Notice that S has the correct asymptotic behavior for 6 — > 0+, with coef- 
ficients which approximate remarkably well the exact coefficients. Notice also 

that the lower bound (Z(^^)(5))" ' is more precise than the upper bound 
Z(^^)(4)/Z(-°^)(5). This behaviour is consistent with Eqs.dM]) and ([gT)) . 

In the case of Neumann boundary condition the exact asymptotic behavior 
for e — >■ 0+ of the fundamental mode of the string (15(11) is not known. It is 
however straightforward to obtain rigorous upper and lower bounds for this 
energy, using the exact sum rules. In Fig. |6]we show the bounds obtained using 
^(^^)(3) and 2'^^^) (4), and the Shanks transformation obtained using using 
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Figure 4: Bounds for the energy of the fundamental mode of the string l|50| l with Dirichlet 
boundary conditions as function of e. 



,(DD) 




4.945 
4.940 
4.935 
4.930 




0.05 



0.10 



0.15 



0.20 



Figure 5: Comparison between Eq.JSTJ (sohd hne) and the Shanks transformation 
[Z(00)(3)-l/32(DD)(5)-l/5 _ 2(i3D)(4)-l/2] / [z(00)(3)-l/3 + ^{DD)(5)-l/5 _ 2Z(DD)(4) 

(dashed hne) as function of e. The lower and upper curves are the lower and upper bounds 
respectively. 
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,(NN) 




4.65 F 



Figure 6: Bounds for the energy of the fundamental mode of the string II50I I with Neumann 
boundary conditions as function of e. 



ZiNN)(^2), Z(^^)(3) and Z(^^)(4). The shaded area corresponds to the allowed 
region. 

If we expand the expression obtained with the Shanks transformation for 
e — !> O"*" we obtain 



{NN) 



4.9336 + e ( 0.7852 cos ( ^) - 0.7852 j 



'27r 

-0.3122 cos ( — 



0.01562 cos ( ^ ) - 0.1084 ) + . . (54) 



Therefore we see that for Neumann boundary conditions the fundamental 
mode of the string is more sensible to the rapid oscillations of the density: as 
a matter of fact, e[^^^ contains a term of order e. This term includes an 
oscillatory contribution cos (^); in the case of Dirichlet boundary conditions 
the dependence on e starts at order e^, while the oscillatory contributions only 
start at order e'^. In other words, it is easier to observe the periodicity of the 
density looking at the Neumann rather than Dirichlet spectrum of the string. 



4. Conclusions 

We have obtained explicit expressions for the sum rules involving the eigen- 
values of strings with arbitrary density for different boundary conditions and we 
have provided simple diagrammatic rules which allow one to obtain the expres- 
sion corresponding to a given order. Despite the factorial growth of the number 
of diagrams to a given order, we have derived general expressions up to order 
9, corresponding to 20160 diagrams, using Mathematica. 
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These sum rules can be used to obtain precise bounds on the lowest eigen- 
value of the string. A more accurate determination of this eigenvalue can then 
be obtained taking into account the known asymptotic behaviour of the spec- 
trum and by performing repeated Shanks transformations of the sequence of 
approximations. Since we deal with exact results, no numerical instability due 
to round-off errors is ever present. In this way we have been able to obtain the 
Dirichlet eigenvalues of a particular string, first discussed by Horgan and Chan, 
with 17 digits of precision. 

For the case of the Borg string, we have proved that the sum rules reduce to 
the analogous sum rules for a homogeneous string only for Dirichlet boundary 
conditions. Therefore the Borg string is isospectral to the homogeneous string 
only in this case. 

For the case of a string with rapidly oscillating density, we have used the 
sum rules to obtain bounds on the lowest eigenvalue: in the case of Dirichlet 
boundary conditions we have verified that the sum rule approximate very well 
the exact asymptotic behaviour of this eigenvalue for arbitrarily rapid oscilla- 
tions of the density, providing the exact functional dependence on the physical 
parameter e. We have then applied the same method to study the lowest eigen- 
value for Neumann be, for which the exact asymptotic result is not available, 
showing that it is more sensible to the oscillations of the density. 

The extension of these results to higher dimensions is treated in a companion 
paper [i]. 

Acknowledgements 

This research was supported by the Sistema Nacional de Investigadores 
(Mexico). 

Appendix A. Green's functions 

In this appendix we derive the explicit expressions for the Green's func- 
tions of the negative laplacian in one dimension and with different boundary 
conditions. 

We need to solve the equation 



with |x| < a/2 and \y\ < a/2. 

Appendix A.l. Dirichlet boundary conditions 

The eigenfunctions and eigenvalues of the negative ID laplacian for Dirichlet 
boundary conditions are 



d^ 

— G{x,y) = S{x - y) 



(A.l) 




(DD) 



(A.2) 



(A.3) 



a- 



,2 
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The Green's function is obtained as 

oo 



ADD) 



(A.4) 



and Ea. (|A.ip follows from the completeness of the basis 
It is easy to see that 



(a + 2a;) (a - 2y) 



4a 



e{y-x) 



(a — 2 max[2:, y]){a + 2 min[a;, y]) 



(A.5) 



To verify this result we just need to check that, for x ^ y, ^ ■^G'"^^^ [x, y) 
and that 



ax 



= 1 



(A.6) 



which follows from integrating Eq. (jA.l[ ) on an arbitrary interval containing y. 
Notice that 



a X 
1~ ~a 



(A.7) 



Appendix A. 2. Neumann boundary conditions 

The eigenfunctions and eigenvalues of the negative laplacian for Neumann 
boundary conditions are 



, n = Q , u — 1 
, n > , u = 1 



INN) 



2 • {2n — l)7ra; \ n o 

- sm -^^ — , n> \) , u — 2 



u=l 



021^2^ , u^2 



(Ai 



The Green's function is obtained as 



ri=0 



(NN) 



(NN) 



and Eq. (|A.l|) follows from the completeness of the basis Notice 
that this expression is formally divergent, because of the zero mode which is 
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present in the Neumann spectrum. However the eigenfunction corresponding to 
n = and m = 1 is a constant. 

The derivation of an exphcit expression for G^'^^^{x,y) requires a careful 
discussion because of the presence of the zero mode. Let us write 



where 



(NN) 



E 



(NN) 



E 



(NN) 



(A.IO) 

(A.11) 
^A.12) 



Now 



dx 



oo 



Six-y)-i.T\x)4r\y) 



Therefore, for x ^ y one must have 

j2 



^^&^^)[x^y) = -1 
dx'^ a 



(A.13) 



(A.14) 



To ensure that the rhs of Eq. ljA.ll) is obtained we also need to impose that 



dx 



dx 



(A.15) 



However since g\^'^^ (x, y) is a constant, the equation above may be cast 
directly as 



It is easy to see that 



(A.16) 



x->3/- 



(a^ + 6a(j/ 


- cc) 4 


6(x2- 






12a 






[a? + &a(x 




6(a;2- 






12a 






{o? ~ 6a|x 








12a 



6{x - y) 
e{y - x) 



and 



a X 
12 ^ T 



(A.17) 
(A.18) 
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Appendix A. 3. Mixed boundary conditions: Dirichlet- Neumann 

The eigenfunctions and eigenvalues of the negative laplacian for Dirichlet- 
Neumann boundary conditions are 



(A.19) 



^i^^H-) = ,n>l 

ADN) ^ (2n-l)V 
4a2 

The Green's function is obtained as 

n=l 

The derivation of the Green's function for this case is completely analogous 
to what done before for Dirichlet boundary conditions and therefore we only 
report the results: 

= (x + a/2)e{y - x) + {y + a/2)9{x - y) 

= {mm[x,y] + a/2) (A.20) 



and 



G(^^)(x,x) = (A.21) 



Appendix A. 4. Mixed boundary conditions: Neumann- Dirichlet 

The eigenfunctions and eigenvalues of the negative laplacian for Neumann- 
Dirichlet boundary conditions are 

= «(^) ,n>l 

e(--) = ^^^^^ (A.22) 
The Green's function is obtained as 

n=l 

Once again we only report the result: 

= {-x + a/2)e{x - y) + {-y + a/2)9{y - x) 

= (- max[a;, y] + a/2) (A.23) 

and 

G(^^)(x,x) = -x+l (A.24) 
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Appendix A. 5. Periodic boundary conditions 

The eigenfunctions and eigenvalues of the negative laplacian for periodic 
boundary conditions are 



n = , u = 1 

^\!^^\x) = { ^/icos^ , n>0,u==l 
^sin22Hi£ n > , M = 2 

\ V " " — 



^n.u — „j (A. 25) 



The Green's function is obtained as 



(PP) Zl^ (PP) 

n— ^nl n— 1 n2 



In this case the same considerations done for the case of Neumann be apply 
and we have 



and 



+ 6a (a; 


- y) + 6(x 


-V? 




12a 




+ 6a(y 


— x) + 6(x 


~y? 




12a 




— Qa\x 


- y| + 6(x ^ 


-y? 




12a 






a 
12 





{x - y) 



(A.26) 



(A.27) 



Appendix B. Asymptotic laws 

We derive here the leading asymptotic behavior of the spectrum of an in- 
homogeneous string subject to any of the boundary conditions discussed in 
this paper. We follow the discussion of Ref. [2^ and consider the operator 

O — whose spectrum coincides with the spectrum of a string 

of density T,{x). 

Let e„ and ipn{x) be the eigenvalues and eigenfunctions of the negative Id 
laplacian: 

3-5- = e„^n{x) ■ 

Define: 



(B.l) 
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where a{x) = J^g^,2 \/^(jj)dy- Notice that 5'„(a;) and ipn{x) obey the same 



boundary conditions only for the Dirichlet case. 
Then 



0*„(a;) 



4I](x)E"(a;) - 5I]'(a;) 



(B.2) 



16S](a;)3 

2 

^(^12)-^ are the eigenvalues of a homogeneous string of length 



where e„ 
cr(a/2). 

Taking into account that In oc r? for n — >■ oo, we see that the 5'n(a;) tend to 
become eigenfunctions of O in this limit; we can therefore introduce the operator 



p = o- 



such that 



and write 



= P 



4S](x)I]"(a;)-5S]'(x)^ 
16I](x)3 ^ 

P^'„(a;) = e„^'„(a;) , 

4E(a;)E"(a;) - 5E'(a;)2 



16E(a;)3 



(B.3) 



(B.4) 



(B.5) 



We can now apply perturbation theory and calculate the corrections to the 
eigenvalues of the string, treating the second term in O as a perturbation: 



= i,n\V\n) 



-a/2 \^Y.{xY 



In the limit n co 
1 



C^(a/2) 7-a/2 



4E(a;)E"(x) - 5I]'(x)^ 



16E(a;)5/2 



da; + O 



We can thus read off the asymptotic coefficients A\ and A2 of the string; the 
calculation of the higher order coefficients can be performed in a similar way. 
For the reader interested in this calculation we refer to Ref . 22| , where we have 
derived the analytical expression for A3 (see Eq.(53)). 

Therefore 



En £n 



1 



a{a/2) J_^/2 



4E(a;)S"(a;) - bY.'{xf 



16E(a;)5/2 



dx + O 



1 



(B.9) 
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